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Free vibration characteristics of a t hin walled 
open cross section beam i.rith unconstrained damping layer is 
studied in tiie present work. A bj/'pical example of a top hat 
section with viscoelastic layer at the flanges is considered. 
Transverse vibrations both in and out of the plane (coupled 
bending and twisting modes) of symmetry are investigated. 
Results are obtained for natural frequencies and modal loss 
factors in both these modes of vibration. Beams -i/ith simply 
supported and clamped end conditions are investigated. 
Extensive parametric investigations are carried out with an 
efficient algorithm based on two dimension Newton-Raphson 
technique . 



CHAPTER 1 


INTRODUCTION 


1 .1 Introduction 

Control of vibration is an important aspect in the 
design of resonant systems. This is necessary in order to 
avoid early fatigue failure and malfunctioning of such sys- 
tems. .Vibrations control is a3.so of significance in view of 
human comfort especially in vehicular systems. 

The level of vibration in a system can be reduced 
in a number of ways. These include, among others, 

(i) reduction of the excitation, (ii) isolation of the source, 
(iii) introduction of vibration neutralizers and (iv) detuning 
and decoupling by proper design. For resonant systems, how- 
ever, increase in the damping is the most effective way of 
controlling excessive vibration. Several damping treatments 
have been devised and practiced towards this end. For 
example, damping tapes are applied to aircraft structures with 
a view to increase the fatigue life and reduce the interior, 
cabin noise. Mastic deadners serve the same purpose in 
automobiles . 

The fuselage of an aircraft is subjected to con- 
vected random pressure field from the jet noise and turbulent 
boundary layer. As a result, -the fuselage vibrates and in 
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turn radiates noise into the passenger cabin. Cabin noise in 
the raid and high frequency regions can be adequately atte- 
nuated by skin damping tape or fiberglass insulation. How- 
ever, these methods are not very effective at low frequencies. 

A concept similar to that of the tuned vibration 
damper has been attempted [ 1 1 in the design of the aircraft 
skin construction to reduce the vibration levels at low 
frequencies. This concept has been referred to as 'intrinsic 
tuning'. For a typical skin-stringer combination, the fre- 
quency of peak response is very close to the natural 
frequency of the individual skin bay (f„)) clamped along the 

Jr 

stringers and simply supported along the frames. If stringer 
and panel dimensions are such that the natural frequency of 
the stringer (f ) concides with f , then the stringer can act 
as a tuned absorber. To increase the effectiveness of the 
intrinsically tuned stringer over a broad frequency range, it 
is highly desirable to introduce additional damping. Addi- 
tional damping in the stringer can be incorporated by various 
methods. The most common methods are 

i) to have unconstrained or constrained visco-elastic 
layers on the stringer flanges, 
ii) to use inherently high damping composites or 
alloys as the stringer material. 

Stringers are thin walled open sections beam type 
structures. No adequate theory is available for the vibra- 
tion analysis of such structures with layered damping 
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treatment. The purpose of the present thesis is to mahe a 
theoretical investigation of such problems. 

1 .2 Short Review of Previous Work 

Work on layered damping treatments was initiated 
by Oberst and Lienard t2] . They considered the problem of 
an unconstrained layer damping treatment. The overall damping 
capacity of the 'composite plate was expressed in terms of the 
loss factor, Kerwin [ 3 ] investigated the problem of a cons- 
trained visco-elastic layer where loss modulus in shear 
contributes predominantly towards overall damping. Ross, 
Kerwin and Ungar [ 4] derived a general theory for free har- 
monic wave motions of an infinitely long beam made of three 
layers . An expression for the overall loss factor as a 
function of frequency was derived. The theory was valid only 
for simply supported end conditions in the case of finite 
beams. Ditaranto [ 5 1 presented a general theory for free 
vibrations of sandwich beams with arbitrary boundary condi- 
tions. Complex natural frequencies ( w ) were obtained in 
the form = w (1 ■f i q ) where n is the loss factor. 

Interestingly, it was found that the modal loss factor 
bears a relationship with which is independent of end 
conditions. Moreover, this relationship was found to be 

the same as that given by the previous work on infinite beam. 

o 

However, the values of w depend on the boundary conditions 



and as such the nnmeriGal values of ri need to he determined 
separabely for each boundary condition. 

Mead and Markus [6 ] considered the forced vibration 
of a three layer damped sandwich beam with arbitrary boundary 
conditions. They used the concept of forced damped modes for 
the response calculation, Markus et al [ 73 also discussed 
numerical techniques for determination of modal loss factors 
and natural frequencies of such structures. The problem of 
damped sandwich plates with elastic constraints along one pair 
of edges was investigated by Mead (81 . He also investigated 
the free vibration of a damped sandwich plate on periodic 
supports [9]* 

All these works, mentioned above, consider beams of 
solid cross section for which the centroid and the shear 
centre coincide. For such structures the bending and tor- 
sional modes are decoupled. 

As already mentioned, stringers are beam type 
structures with open cross sections for which the shear centre 
and the centroid usually do not coincide. In such cases 
there will be coupling between the bending and torsional 
modes of oscillation. Gere and Lin [10] presented an approxi- 
mate method based on Hayleigh - Rltz formulation for deter- 
mination of the coupled natural frequencies. 
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The major objective of the present work is to 
provide a theoretical analysis for the free vibration charac- 
teristics of thin walled open section beams with -unconstrained 
viscoelastic layer. As an example, the problem of a top hat 
section with damping layer on the flanges is considered. The 
base material is ass-umed to be nondissipative and the loss- 
factor of the viscoelastic material is assumed to be indepen- 
dent of frequency. Moreover the loss factor is assumed to be 
the same in the extensional and shear modes of deformation. 

The exact analysis of the free vibration characteristics of 
thin walled open sections is presented in Chapter 2, Numeri- 
cal res-ults of natural frequencies are co-mpared with those 
obtained by the approximate method of Gere and Lin [101. 

The analysis for the uncoupled bending vibration in 
the plane of symmetry is essentially that of the transverse 
vibration of beams with solid cross-section 12 to 9 3* The 
natural frequencies and modal loss factors are obtained in 
close form in terms of geometric and material properties of 
the beam. 

The analysis for the coupled bending and torsional 
oscillations require, as a first step, the determination of 
the St, Venant's torsion constant and the warping rigidity 
of the composite beam. These quantities turn out to be 
complex quantities due to the presence of the dissipative 
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layer. The St, Venant torsion constant is determined by 
extending the membrane analogy of elastic sections to the 
composite secbion. The warping rigidity is derived in terms 
of the average warping function. 

The coupled bending and torsional modes are gover- 
ned by coupled eighth order differential equation. Assuming 
harmonic motion for the free vibration, the frequencies turn 
out to be complex. The real part is the natural frequency 
and the ratio of the imaginary part to the real part gives 
the modal loss factor. 

The standard procedure of determining the frequency 
results in a frequency determinant with complex elements. 

The zeroes of this determinant determine the natural frequen- 
cies and modal loss factors. The complex arithmetic adds a 
dimension of difficulty to the determination of roots of the 
frequency equation. This difficulty is overcome by using a 
two dimensional Newton - Raphson technique. Numerical 
results are presented for various geometric and material 
parameters involved. 



CHAPTER 2 


FREE VIBRATION OF OPEN SECTION THIN WALLED BEAMS 


2.1 Formulation of the Problem 


Figure 2.1 shows a typical beam cross section 
considered in this thesis; where 0 is the shear centre and 
C is the centroid. X axis is the longitudinal axis through 
0 forming a right handed coordinate system with Y and Z 
axes. In the following analysis v and w represent 
deflection along Y and Z axes respectively, refers to 
the rotation of the cross section about X axis . The coor- 
dinate system xnC iias its origin at C. 

It can be readily seen that 


where c- is the distance of the centroid from 
z 

the shear centre, the subscripts refer to the points under 
consideration. 


E 

E I^ 


The eauation of motion can now be written as 1 101 


h .2- 

g Wq g 

-r— + P — p 

gt^ 

+ p A 

g X 9 1 


ax 

g V 


= 0 


( 2 . 2 ) 


2 (Vq - ) =0 


+ P A 


(2.3) 




9 


E C 


3 ij; 


w 


3x 


If" 


- G C 


® ax 

+ Pi 


p A ( 
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(2 .If) 


where and are area moments of inertia of the 

cross section about ri and ? axes respectively* A, C_ and 
are area of cross-section, Saint Venant's torsion constant 
and the warping constant. E and G are elastic and shear 
moduli of beam material. is the polar moment of inertia 
of beam cross-section. 


Equation (2.2) is the classical Euler's equation 
for vibration of the beam for which the results are well 
known . 

Assuming harmonic motions for the free vibration 


V , = V (x ) e^ “ 

o o 

and ip = 

equations (2.3) and (2.!+) can be transformed, using equation 
(2,1 ) to 
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( 2 . 5 ) 
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In terms of non dimensional quantities, equation 
(2.5) can be rewritten as 


— -o + a-1 YZ + 


2 

d^ V 


ou ■§-+n'^a, V 

n3^^2 n4o 


( 2 . 6 ) 


The constants etc, are given in Appendix 

A, and is the nondimensional natural frequency, 

8 \ X 

Assuming the solution V (x) = V e 

for equation (2,6), the follov/ing characteristic equation is 
ob bained for x 


8 


+ 


+ Ct, 


Q 


1 x'* 
n 


+ a. 


n 


+ a 


Q 


4 


\ . it 

4 n 


= 0 


(2.7) 


Applying appropriate boundary conditions at the 
two ends of the beam the frequency determinant (order 8x8) 
in the form 

Det t A] = 0 (2.8) 

is obtained. 

The elements of the matrix A for simply supported 
and fixed end beams are given in Appendix A . The coupled 
natural frequencies, ^2 n are obtained by using Newton - 
Raphson method. 


2.2 Results and Discussion 

Table 2,1 shows the first eight natural frequen- 
cies for the coupled vibration. The table also includes 
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the results obtained by the approximate method of Gere and 
Lin [10 ] , wliich are found to be in close agreement. Table 
2.2 presents the natural frequencies in the vertical (uncou- 
pled) and in the coupled modes, nondimensionalised with 
respect to the same quantities. It is seen that some of the 
natural frequencies in the vertical mode are quite close to 
those in the coupled mode. This suggests the possibility of 
these modes being excited simultaneously. This behaviour has 
also been observed experimentally [11 ]. 
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TABLE 2.1 


COUPLED NATURAL FREQUENCIES BY EXACT AND APPROXIMATE 

ANALYSIS 


= 2.^00 cm 

= 1 .7780 cm 

h = 2.9718 cm 

t = 0.2032 cm 


Q 


n 


_ “n 
" **1 1 

I 

P A 


Q 

a 



where a3„ and w are the natural 
n a 

frequencies by present analysis and 
approximate analysis respectively. 


END CONDITIONS 


MODE 

1 Simply Supported Ends 


1 Fixed Ends 

1 Natural Frequency by 

1 Present Analysis 

5 

5 Natural Frequency j! Natural Frequency 
i by Approximate by Present 

I Analysis [ 10]| Analysis SI 

1 

7.40 

7.40 

10.13 

2 

17.92 

20.86 

26.81 

3 

20,87 

23.91 

46.23 

k 

38,71 

47.79 

84.70 

5 

67.21 

74.76 

126.98 

6 

81 .60 

84.92 

174.91 

7 

104.23 

109.67 

232.42 

8 

149.40 

174.12 

248.54 
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TABLE 2 ,2 

NATURAL FREQUENCIES IN COUPLED AND VERTICAL MODES 
OF A SIMPLY SUPPORTED BEAM NONDIldEl'ISIONALISED IN 
H-IE SAME MAI^INER 


2 . 5^00 

cm 
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0 

II 
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a - - 




1 .7780 

cm 
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PA 1 

2.9718 


cm 





where 

and 
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0,2032 

cm 

c 
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frequencies in coupled and ver- 
tical. modes respectively. 


MODE 

5 Frequency in Coupled 

5 Mode, 

5 ° 

Q Frequency in Vertical 
I Mode , 

5 V 

1 

8,61 


9.88 

2 

28.60 


39.51 

3 

33.27 


88.89 

1 + 

61 J+6 


1 58.0if 


107.26 


246.93 

6 

130,23 


355.58 

7 

166.3? 


483.99 

8 

238 »4-4- 


632.16 



CHAPTER 3 


EFFECT OF UNCONSTRAINED DAMPING LAYER ON . 

GEOMTRIC AND ELASTIC PROPERTIES 

3»1 Shear Centre Location 

When the damping laj'^er is a b bached to the flanges 
of the cross-section as shown in Fig. 3*1 the shear centre 
of the composite cross-section will be shifted. The new 
position of the shear cenbre is determined by bhe following 
analysis . 

Consider the bending of bhe section about the Z 
axis under the application of a shear load V in the Y direc- 
tion through the shear centre 0. AA. and BB are two adjacent 
sections along the length of the beam, distance dx apart , 

1 and 2 refer to the base material and the viscoelastic 
layer respectively. Consider a section CC at a distance y 
from the neutral axis Z . The unbalanced axial force due 
to the differential bending moment is eq'uilibrated by the 
elemental shear force dF given by 

^0 

dF = 1^ J ( dz 

^o 

z^+t^ +t2 

, + 1-j J ( "■ (3.1 ) 

where a refers to the bending stresses, with the subscripts 
indicating the material and section. 




Fig. 3-1 Section of the composite beam 
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The bending moment at a section is resisted by both 
the base layer and the viscoelastic layer. Assuming the 
longitudinal strain at a section distant y from the neutral 
axis is the same for both, the portions of the bending moment 
resisted by each can be shown to be 

M E. Ir 
‘ 1 

M = (3.2) 

H EJ 

^2 = ^ — i 0.3) 

where M is the applied bending moment and j ? 
the area moments of inertia about the neutral axis, and E^ and 
E 2 are the elastic modulus of the base material and the complex 
modulus of elasticity of the viscoelastic material. 


FromEqns. (3.1 ) to (3.3) and using engineering 
bending theory it can be shown that the shear stresses in 
the two layers are 


V E^ 

V E* Qp 


(3.^) 


where and are the static moments about the neutral 
axis of the areas above the section at which the shear 


stresses are acting 
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Considering the moment balance about a convenient 
point in the plane of the cross-section of all the shear 
forces the shear centre can be located. The final expression 
for the shear cenbre location is given in Appendix B . It 
should be noted that the shear centre position in this case 
is not only a function of the geometrical properties of the 
section but also depends on the material properties. Because 
of the complex modiiLus of the viscoelastic layer the location 
is given by a complex number. 

3*2 Torsional Rigidity of the Two Layer Beam 

In the case of thin walled open section beams, the 
applied torque is resisted by two different types of defor- 
mation. The first part corresponds to the torque balanced 
by shearing stresses due to pure torsion (St. Venant's torsion). 
The other part is the torque resisted by the warping stresses.- 

The St. Venant's torsional rigidity corresponding 
to the torque resisted due to pure torsion for the composite 
section is derived by an extension of the membrane analogy 
method. 

Consider the equivalent membrane problem, namely 
that of the deflection of the membrane made of two materials 
of equivalent size (Big. 3«2a) subjected to internal pressure. 
The membrane w lich is initially in the form of 0 = 0 will 
deform in the shape 0 = 0 (z, y) as shown in Figf 3»2 b. 
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Fig. 3-2a Equivalent 
membrane 


Fig.3-2c Approximate shape of 
deflected membrane 


Y 

/ 



Fig. 3-2 b Actual shape of deflected 
membrane 
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If 1^ >> and >> t^, the cross section 
of the membrane can be considered as independent of y 
except near the end as shown in Fig. 3 •2c. 0^ and 02 are 

the deflected surfaces of the membrane of material 1 and 2 
respectively which also correspond to the Saint Venant's 
stress functions. 


Let 0^ be the unifom twist along the length, 
then using the membrane analogy 0^ and 02 are governed by 
the following equation, 


8^ A 


3 z 
3^ 0 


■* 2 G, 


(3.^) 


a z 


r = - £ G2 *0 


where and G 2 are respective shear moduli 


^’’'xy 


The shear stresses can be determined from the stress 

9 00 


functions using the relations 

a 0 




az 


and 




a z 


(3.6) 


0^ and 02 can be found by integrating Eq. (3*5’) 
and applying the following boundary conditions s 


0^ (0) = 02 (0) 

0^ (- t^ ) = 02 ('t2^ - G 


At 


i.G. 


= 0 , 


'V’l 


'^xy^2 


a 0 


1 


1 


1 


3 z 


z=0 


Qf 




3 02 

a z 


z=0 


where 


Y represents the shear strain* 


(3.7) 
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^1 

^2 


The final forms of and ^2 turn out to be 

= “ 0^ (2^ - tp + (z + t^ ) 


(3.8) 


” ^2 ” ^2^ '*■ ^^2 ®o - tg) 


where U 


G* t| - tf 


to + G^ t. 


2 2 


'1 "1 


The torque resisted by the shear stresses in the 
flanges induced by pure torsion is given by 


T^ = 2 I J 0 dz dy 

= 2 / d + 2 / ^2 ^ ^2 

^1f ^2 

where and A2 are tho areas of cross-section 

of the flanges for material 1 and 2 respectively. The 
torsional rigidity of the entire composite beam including 
the torque resisted by the rest of the cross-section is 
given in Appendix B 


3*3 Warping Rigidity of the Two Layer Beam 

The warping stress generated due to fibers non- 
uniform torsion, induce shearing stresses which results in a 
re sis fcive torque . 

The warping displacement w^ can be written as [I7] 


w- 


’'s - 


''f 


.aiiL. 

9 X 


(3.10) 



where 


Average warping function, 
Warping function. 

Rotation of cross-section. 
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w 


w. 




Since ^ varies along the length of the h 


earn, 


adjacent cross sections will not be warped equally and there 
will be axial strain e in the longitudinal fibres of the 

ji. 

beam, given by 
3Wj 


'X 


3x 


( - w„ ) 


d^l p 

3X^ 


(3.11 ) 


Assuming no slipping at the interface, the axial 
strain should be the same in both elastic and visco-elastic 
parts. But the warping stresses produced due to the axial 
strain are different due to different elastic inoduli of the 
two materials . 


Assuming there is no lateral pressure between the 
longitudinal fibres, one can write 


CT 

X 


1 


) 


d^lp 

3x2 


and 




( w 

2 S 


W^ ) 

s 


3x2 


( 3 . 12 ) 


where a. 


X. 


and 


X, 


are the warping stresses in the 


two different layers of elastic moduli E-j and E 2 ' 


Figure 3.3a shows an element cut out from the 
flange of the elastic portion of the beam shown in Fig. 3.3b. 
Along the middle line, the stresses vary with distance s 
measured from the edge of the section and can be calculated 
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from an equation of static equilibrium of the element 

m n o p. Projecting all forces on to the x - axis, one can 
wribe, 

9 ( 'T-] t. ) 

^ — ■ ds dx + t^ — g— ds dx = 0 (3.13) 


Integrating with respect to s and observing that 
^ is independent of s and that vanishes for s = 0, 

the expression for can be wihtcen by using equation (3-12) 




9 x-^ o 


1 ( 


Wg - Wg ) t^ ds 


( 3 . 1 ^) 


The resistive torque is obtained by summation along 
the middle line of the section of the moments of the elemental 
shear force about the shear centre. Therefore, the expression 
for the warping rigidity for the elastic portion of the 
beam can be written, using equation (3.1^) 


M 


1 



9^ 


9 X' 


ra s 

i ' i 


( w - w ) t. ds } r ds 

O O 1 

( 3 . 15 ) 


where m is the perimeter of cross-section and r is 
the perpendicular moment arm about the shear centre . 


equation 

equation 


Integrating the last part of the integral of 

m _ 

(3.15) by parts and using J ( -w - w ) ds = 0 , 

o s s 

(3.15) can be written as 

9^4; _ Q 

= - E. o / ( w - w ) t. ds 

‘ 9x-^ o ® 


M 


-^1 ^w 


934; 


1 ax 


3 


( 3 . 16 ) 
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where C,, is the warping constant for the elastic 

W w| 

portion of the heam as given in Appendix B, 

Similarly, one can proceed to find the expression 
for Mg, the warping rigidity for the visco-elastic part, as 


14 . 


- Eo C 

2 w. 
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3x 


3 


(3.17) 


where C is oho warping constant for visco-elastic 


w 


portion of the beam and is given in Appendix B. 


Table 3*1 presents some calculated values of the 
quantities discussed in the Chapter for a Typical Cross- 
Section . 
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TABLE 3.1 

SHEAE CENTRE LOCATION, TORSIONAL RIGIDITY AND 
VJARPING RIGIDITY V^ITH VISCOELASTIC LAYER 


— = 0,001 , 3 = 1.0 

h = 3.50 cm, 

E^ = 70 X 10 ^ Kgf/cm^ 


=3*0 cm, I 2 = 2,60 cm_, 

= 0.20 cra^ 


t5 5 Shear Centre Location ( Torsional Rigidity ' "II Warping Rigidity 

5 ^le ) 5 , Cl , 5 . . C^ 

1 iSReal Part S imag. Pa'rtjReal Part ft Imag. Parti Real Part pmag. Part 
I cm I cm ^ ^ in-“2 5 1 n-^ 5 

I jj I Kgf -cm I Kgf - cm jKgf - cm^ | Kgf -cm^ 


0.50 

0.71677 

-0.00076 

0.10165 

0.10443 

0.70857 

0.17797 

1 .0 

0.71599 

-0.001 5if 

0.101 93 

0.38755 

0.70876 

0.36402 

1 .50 

0.71520 

-0.00233 

0.10249 

0.93828 

0.70896 

0.55827 

2-0 

0 . 711 + 4-0 

-0.00313 

0,10340 

1 -8453 

0.70916 

0.76085 



CHAPTER If 


FREE VIBRATION OF OPEN SECTION THIN WALLED BEA14 
WITH UNCONSTRAINED DAJIPING TREATMENT 
AT KIE FLANGES 


If.1 Free Vibration in the Plane of Synmietry 
If, 1.1 Eqiiation of Motion 

Figure If,1 shows the beam cross section with a 
visco-elastic layer on the flanges is the neutral 

axis of the untreated beam section, is that of the 

visco-elastic layer and n - n is the composite neutral axis. 


For displacement in Z direction the strain at 
and ^2 " expressed as [^3 


^ D 


8^ w 


a X 


H, 


8^ w 


(If.1 ) 


20 3^2 


where x is measured along the length of the 
beam and other quantities are explained in Figure 4,1, 


The net extensional force on each layer can be 
obtained from tlie strains at their respective neutral axes 
as 




^2 "^2 ^^2 0 


8 X 


and 


(4.2) 
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As the total extensional force on the composite 


section is zero, one gets 


+F2 = (E* AjHgo - E, D) 

9 X 


0 


or 


^2 ®21 “ D 


Thus the shift of the composite neutral axis from 
that of the original beam cross section, D, is given by 

p = + i -Dj (say) ... C^.3) 

^2 A^ + E^ A^ ^ 

with E 2 = E^ (1 + i 3 ), where 3 is the loss factor of 
the visco-elastic material, one gets 



E A 
_2 ^2 


(E^ A^ + E 2 A^)^ + 3^ e| a| 


(4.4) 


and 


D-, 


3H2^ E^ A^ 1&2 ^2 


(E^ a^ + E 2 A2)^ + 3^ e| a; 


( 4 . 5 ) 


The equation of motion can be derived by using the 


relation 

ill = - (P. A. + Pp Ap) -i-W (4.6) 

3x2 114 6 5^.6 

The expression for the bending moment M can 

be written as 


M 


= *Sl +«22 +^’2H20 
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vhere M, 


•^2 

1? T d W 

•^1 -^n p 

'l1 9 


and M, 


s2 

d w 


F* T -11- 

•^p -*-n 2 

^ '22 9 


(l-F.?) 


with I and as the second moments area 

^1 1 *22 

of the layers about their own neutral planes. Using equa- 
tions (4.2), (^.6) and (4,7) the final equation of motions 


is obtained as 


B (1 + i V ) 

9 x^ 


- p A -'f I 

9t 


(4.8) 


where p A = P-j -4^ + A 2 


B (Df-Df)} I + 


Aj ( h|^ + d| - d| - 2 Dj, H 21 - 2$ Dj (Dj, - > 

(4.9) 

and 

V B = 2 A^ Dj + 2 A^ E2 Dj 

6^2 < ^22 •^'^2 ®ll 

(4.10) 

Now, equation (4.8) is simply that of a hyst^ritically 
damped beam. So natural frequencies and loss factors at 
different modes for various boundary conditions can simply 
be written as follows: 
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For simply supported ends 


w 


If If 
n TT 


B 


n 


1 


T" pA 


n = 1, 2, 3, 


and p 


For claraped clamped ends 
, (n + lA 


(Jj 


n 


1 


IT 


B 

P A 


and n = V 


(^,11 ) 


(If. 12) 


If Eg << 5 as is usually the case, B and v 

can he approximated from Eqs, (4.9) and (4.10) as 




(4.13) 


E 


+ A. H 


and n = V 3 


_ 2 . 
E 


22 


2 21 


1 ^ 


(I+.14) 


4.1 .2 Numerical Results and Discussion 

It can easily be seen from Eq, (4.11) and (4.12) 
that tlie loss factor in this mode of vibration in the plane 
of symmetry, is independent of both the end conditions and 
the mode number. 

Figures '4.2 and 4.3 show typical values of the 
loss factor for various values of tg/t^, Eg/^-] 
loss factor of visco-elastic material 3 . The following 
dimensions are used for the cross section of the elastic 



Loss factor (f\ } 
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V n- 0'6 
0 p 1-0 

• n 20 


E2_ 


0-001 


P2_h 



Thickness ratio (t 2 /t^ ) 

111 } 

Fig.4‘2 Variation of loss factor with 

vertical mode of vibration “““ __™ 

Hitt. ■MS 
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layer, 1^ =2.00 cm, = 2.60 cm^ h = 3*5 cm and 

= 0,20 cm. As expected from the approximate expression 
(^-.1*+), the composite loss factor n is seexi to be almost 
proportional to and also proportional to 3 . The 

nature of variation of ri with t2/t^ is also as predicted by 
Eq. 0-h.^h-). 

4-. 2 Free Vibration in Coupled Bending Torsion Mode 

’+.2,1 Equation of Motion 

Consider the top hat section with the damping 

layer in its flanges as shown in Figure 4.4-. 0 and G 

respectively are the shear centre and the mass centre of 

composite section. Let v and v denote the displace- 

V w 

raonts of C and 0 along the Y axis respectively. It 
is obvious that 

Vc = Vo - C^ if) ... (4.15) 

where is the distance of the mass centre from the 

shear centre. Considering the equilibriimi of elastic and 
inertia forces in Y direction in one hand and the moment 
equilibrium about a longitudinal axis through the shear 
centre on the other, the following equations of motion 
are obtained. 

4 2 
* S ^ 8 V 

V + 4 t ^ “1 4 + Oz 4^ --r = ° 

2 9t 

(4,16) 




Fig.4-4 Cross seetion of composite beam 
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'2 


G, 




ax 


2 ( p-j A-| -f p2 A2 ) 


i^Zc c 
at^ ^ 


+ ( 0^ I,, + P 2 le ) 


= 0 


(^.17) 


■where and are Young's modulus of elasticity, 

and are area moments of inertia about t axis, 

2 ^ 
and Ag are areas of cross-section, and are mass 

densities, I and I are polar moments of inertia about 

a longitudinal axis. Subscripts 1, and 2 refer to the base 

material and viscoelastic layer. , C2 are the complex 

torsional rigidity and warping rigidity respectively of the 

composite cross section as derived in Chapter 3* 


Assuming harmonic solution of the form 


V = V (x) 

o o 

^ (x) 

and using the relation for eliminating ijj , the equa- 

tion of motion can be transformed into the eight order 
ordinary differential equation in nondimensional form as 


d® V 


d 


o 


‘1*’ '^o 2 A 

+ a. - — a 




+ 


c 3 d 


+ 



a I, 



(4,18) 


where ct^, a^, c^, and 2 ^ are given in Appendix C. 
The final form of equation (4.18) is exactly similar to 
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the equation (2.6) for the beam without the layer- But 

“l5 ^ 2 ’ ^3’ complex quantities in this 

case due to the presence of viscoelastic layer. In order 

to avoid confusion, let these complex natural frequencies be 
>i:p P where ^ 

denoted as = (1 + ir)) /the real part is referred 

to as the natural frequency and n is the modal loss factor. 

4.2.2 Frequency and Loss Factor Determination 

Assuming a solution of the form 
8 X _x 

Vq (x) = I Q e ^ 

^ s=1 ® 

and substituting into the differential equation the follow- 
ing characteristic polynomial equation for ^ is obtained 

. ^2 .4- 

X^ + ct, X^ + ^ ^ ° 

(4.19) 

Applying appropriate boundary conditions at the 
two ends of the beam the frequency equation is obtained in 
the form of 

Det [A] = 0 (4.20) 

where the elements of matrix A are given in 
Appendix C. 

The solution of the determinental equation in ' 
this case is much mote complicated as compared to that 
without the viscoelastic layer • This is so because of 



the cumbersome complex arithmetic involved and the inherent 
numerical problems associated vith it. This difficulty is 
overcome by splitting the determinant into its real and 
imaginary parts and setting these separately equal to zero. 
The real and imaginarj/" parts of the determinant are consi- 
dered as function of tuo variables , n namely the 
real and imaginary parts of the complex natural frequency 
squared. The zeros of the functions are obtained by employ- 
ing a two dimensional Newton - Raphson procedure . The itera- 
tions in the Newton - Raphson procedure proceeds according 
to the following sequence. 

Let 

(x, y) = 0 (4.21) 

Aj (x, y) = 0 (4.22) 

2 2 
where x = arid, y 


Aj^ and Aj are the real and imaginary parts of 
determinental equation. Then, the iterations for the roots 
of Equations (4.21) and (4.22) are 


(x) . 


3+1 




"* hi h. 


- (r) + ( ^ ^ — ) 

3 A„ A.^ - Apj^ Aj^ j 




'X 


and (y )j_j,'j " 

where suffixes 
with respect to 


(y)j + 


h h 


X 


Aj A« 




It. 


h 4 A 


X and y indicate partial differentiation 
X and y respectively (for example 



_ ^ 3Ajj 

“ 3 X ’ ~ on). The suhiscript 3 

X y 

indicates that the concerned functions are evaluated at the 
point (x, '-Hie required partial derivatives in the 

iterations are computed by a central difference scheme. 
Starting from the lower end of frequency scale the itera- 
tions are continued to obtain the natural frequencies of 
various orders successively. The convergence criterion for 
each natural frequency and loss factor are specified in 
terms of the percentage change in the values in consecutive 
iteration and the absolute value of the determinant. 

Tlie above scheme converged in most of the cases 
within six or seven iterations. The initial starting point 
for the variables did not present any problem in almost all 
the cases. However the determinant values vary over a wide 
range, with frequency. This presented some numerical prob- 
lems of dealing with very large numbers . This problem is 
alleviated by suitably scaling the determinant value at 
tho begining of the iteration scheme for the determination 
of a particular* modal frequency. 

4.2.3 Numerical Results and Discussion 

Numerical results are obtained, using the analysis 
presented in 4.2.2,, for the following dimensions and pro- 
perties of the elastic layer: 

1.^ == 2.00 cm, Ig = 2.60 cm, t,^ = 0,20 cm and 

h = 3.50 CHij P,, = 2.7x1 0'^ i-Cgf - sec Vem^ 

= 70 X 10^ Kgf /cm^, G., = 28 x 1 0 ^ Kgf/cm^, 

= 62.00 cm. 


1 
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Table 4.1 presents the first eight natural 
frequencies and nodal loss factors for a set of typical 
values of the parameters involved. It is seen, unlike in 
the case of vibration in the plane of symmetry, that the 
loss factor depends on the end conditions as well as modal 
number. This is due to the shearing of the viscoelastic 
layer caused by the coupling of bending and twisting modes. 
This fact has also been observed in the case of beams with 
solid cross-sections having constrained damping layer. 


For both types of end conditions the loss factor 
in the first mode is seen to be considerably higher than 
those at higher modes. No definite trend in the variation of 
loss factor is observed with the modal numbers. At this 
stage, no attempt has been made to have physical justifica- 
tion for this behaviour. The mode shapes are very complica- 
ted to visualize any pattern of deformation for the damping 
layer. Hence it is suggested to calculate individually the 
loss factor in the required mode. 

Figures 4.5 and 4.6 show the variation of loss 

tp 

factor with thickness ratio (^) for typical ratios of 
Ep Pp 

^ and and B for simply supported and clamped ends. 

It is seen that in the case of coupled modes also the nature 
of variation with is very similar to that observed in 
the vertical mode (shown in Fig. 4.2). The modal loss 
factor also appears to vary linearly with 3 , the loss 

factor of the viscoelastic material. As already obseryed 
in Table 4.1, the loss factor in the first mode always 



TABLE 4,1 


NATURAL FREQUENCIES Mm MODAL LOSS FACTORS IN 
COUPLED VIBRATION 


= 1 . 50 , 

Pp 

"1 

0.20 

^1 

= 2,00 cm, 

I2 = 2.60 cm 

= 0 . 001 , 


1 .00 

^1 

= 0,20 cm, 

h = 3 - 5'0 cm 


The nondimensional coupled nabural frequency 
2 ^ 

SI ~ — where w is the natural frequency 

( A^ ^ 2 ^2 ^ ^ 

n = the composite loss factor. 


5 SIMPLY 

MODE I 

NO i U 

5 _c 


SUPPORTED 

5 

L... „ H 


ENDS I 

I ^ 

I f 


c 


FIXED 



ENDS 

n 


1 

3.87778 

0.00350 

7.76462 

0,001 72 

2 

13.37329 

0.00175 

20.23136 

0.00129 

3 

16.86532 

0,00124 

39.20238 

0.00108 

4 

29.11707 

0.00126 

65.61231 

0.001 03 

5 

51 .14014 

0.001 06 

103.62929 

0.00098 

6 

66.50152 

0.001 01 

133.62726 

0.00090 

7 

79.^7019 

0.00097 

177.67156 

0.00088 

8 

114.04927 

0.00092 

202.95039 

0.00096 



burns out to "be hir^lier than that in the second mode. Coni" 
paring the valuos of loss factors from Figs. 4.5 and 4,2, 
(for the samo values of the parameters involved), it can 
be said that the unconstrained layer is more effecting in 
damping the lower order coupled modes. 

From Figs. 4.7 and 4.8 it is seen again that 

Ep 

the modal loss factor varies linearly with . This 

trend was also observed for the vibration in the plane of 
symmetry (shown in Fig. 4.3). 

Figures 4.9 and 4.10 show that the loss factor 

Pp 

is virtually independent of — in its feasible range. The 

Pi 

obvious suggestion is that one should go for light visco" 
elastic materials so that the increase in weight is minimum 
with the damping layer. 

Figures 4.11 to 4.16 show the changes in the 
natural frequencies with various parameters. 

It can be noted that the natural frequency is 
virtually independent of 3 . For reasonable values of 
( « .001), the frequency is lowered with higher thickness 

ratio as the mass is increased without almost any change in 

Ep 

bending rigidity. However, \d.th ^ ~ The increase 

in mass can be offset by the increase in rigidity so that 
the natural frequency increases with increasing thiclmess 
ratio. Figures 4.15 and 4.16 suggest the obvious results 
of lowering natural frequencies with increasing thickness 
ratio, the effect being more pronounced with higher values 
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Fig,4-12 Variation 'Of natural frequency with thickness ratio 
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Fig.4‘14 Variation of natural frequency with thickness ratio 
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CHAPTER 5 


CONCLUSIONS 

An analysis for free vibration characteristics of 
a tain walled open section beam vrith unconstrained visco- 
elasbic layer at the flanges is presented. Transverse 
vibrations both in and out of tne plane of syTiiraetry are 
considered. In the later case, the bending and the twisting 
modes are coupled. The natural frequency and modal loss 
factors, for both type of •'/ibrations, are reported for 
various values of the parameters involved. 

The nature of variation of the modal loss factor 
with different parameters are found to be similar in both 
cases. The important observations inc3-Ude the almost 
linear variation of the modal loss factor mth the ratio of 
elastic moduli and with the loss factor of the viscoelastic 
material. However, in the coupled modes, the loss factor 
even with the unconstrained layer, depends on the end 
conditions and the modal number. The treatment appears to 
be more effective in damping out the lower order coupled 


modes . 
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APPEllDIX A 


DETAILS FOR THIN WALLED OPEN SECTION 
A-1 Secbion Properties 


Various dimensions of the cross section are shown 
in Figure A-1 where 0 is the shear centre and C is the 
centroid. 


(i ) Area A ~ 2 t^ + ig t-j + 2 h t^ 

(ii) Centroid position, z = ( (Ij, t^ + 2 h^ t^ 


+ 2 1^ t^ (2 h + t^ ) }/ 2A 


(iii) Area moment of inertia about n axis. 


I, 


1 

ig t^ 


+ I t^ I t^ h (h + ) 


■p " 3 

- a1 

(iv) Area moment of inertia about ^ axis 


I, 


, 1 . n , Vti (1^ 


12 


+ ;^ t^ 1::| + 2 


■1 


^ ■*^1 ^2 
+ ~r~ ^2 


(v) Saint Venant's torsion constant 
C : 


2-1 ^3 . 2 -u +.3 , 1. -1 4-3 


3 ll tf + I h tf 3 I 2 t-^ 


(vi) Location of shear centre 


t. 1 


h^ 

( + ---i-g-- -f ) 

2 


1, 1^ h if b 


1 -"2 ^ 


2 t^ h 1-1 ^ Jh. ) 


“¥ ^ “3 
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(viii) Warping constant 


yhere m ^ai^+l^+ah 

,. , l/-\ h. 0 ~l „ 

~ ~2 l-l (h. + e) g- — 

I. = (h + e / 3 


_ 2 

m w ^ 
s 


^2 " T2^ (h + e)^ + 1^ 

- 6 h. I 2 (h + e) } 

1 3 ~ ^P2 'f 3 O’ + 4 e I2 •" 6 I2 6a} 


whore a = Ig h - 2 1^ (h + e) 

-l(. " ^2 ^ ^ + 3 h Ig h + Ig > 

where b ~ Ig h - 2 (h 4- e) - 2 e Ig 

= -^{ 3 + (h + e)^ 1^ + 3 (h + e) d } 

where d = 1^ (h + e) + e I 2 ^ ^2 


A-2 Won Dimensional Terms in Equation (2,6) 


(1 + + s^' ) 
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S*l 1 


and a 
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"W 
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°z 
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C ^ A 
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1 ia the length of the beam. 
Non dimensional natural frequency 




2 _ 
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w 


E I 




p A 1 

where w is bhe natural frequency 


A-^ 


A. 3 Appropriate Boundary Conditions for the Coupled Vibration 


For simply supported ends the boundary conditions 


are 

V =0 

o 

li 

V = 0 

o 

^ - 0 or 

ii 

and t); 0 or 

For fixed - fixed ends, 

= 0 

= 0 

= 0 or 

=: 0 or 


V. 


0 


V 


o 






at X 
and t 


= 0 
o 



the boundary conditions are 


0 

1 


= 0 
o 

= 0 
O 


at X = 0 

and X = 1 



A-5 

A.^ Elements of Matrix A in Equation (2.8) 


Elements of the matrix A is different for diffe- 
rent boundary conditions. For simply supported endss 
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1 .0 

1 .0 

1.0 

... i.o' 


e 

e'2 


... e8 


> 2 

x2 

^2 

X2 

3 

... X2 

[ A] = 

, 2 

X ^ e 


X2 e^3 

3 

- 


^ 1 

^ o 

2 

x^ ... 

3 

... 


^>-1- 
X ^ e 

^2 ® 

x^t ... 

3 

••• 


x6 

X^ 

'^2 

x6 

^3 

... x6 


,6 

X^ e 

6 ^2 

Xg e ^ ... 

... x6 


Here, ?i| , Xt, ... Xg are eight roots of the characteristic 
polynomial equation (2.7), 



For Fixed ends beam; 


X-| e 


e 


X^ e 


Xg e 


L i|. 
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x’^ a ^8 
Xg e 




Xq e 
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X? » 8 
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APPENDIX B 


geoketeig and elastic properties for composite 

GROSS SECTION 
P .1 Secbion Properties 


Fi.cfure 0,1 shows the composite cross section 
wlioro and " ^2 are the centroidal axis of parts 

1 and 2 res 2 :)ectively, C is the mass centre of the 
compos! be sec bion, 

T - ^ ^ ^ 2 .. ^2 \ 

“ ”^1 *^1 -^2 ^2 

where are given as A and X respectively in 

Appendix A. and P 2 are the mass densities of the 
material of 1 a,nd 2 respectively. 

t2 

A 2 ~ 2 1^ ^2^ ^2 ~ + t^ + 2 ^ 


^1 = 



where 
of part 1 




given in Appendix A. 

3_3 + ( 1 ^ + 1 ^) 

and Ir are the area moments of inertia 
^2 

2 ahont the axis of symmetry 


B .2 Shear Centre Location for Composite Beam 
e ~ "^3 ^2 ” "^1 *^2 h + t ^ + tg ) ~ 2 '?2 ^ 




whore 


B»3 


V, 




(E, I 


■\ \ + En ) 

‘ ? -1 2 5 " 


V. 


and V 


1 


witli E, 


* 


2 1 ? ^1 

(E^ I + Ep I ) Tj: + -3 ) 

1 ^ ^2 

^2 ^2 2 Ip 1 . 

Ep^ (1 + i 6) 5 g- xoss factor of the 

viscoelastic material. 


B»3 l^orsional Rigidity of the Composite Beam 


Cr^ ( I t^ + I U t^) 

+ h G* { I J y t|) 

+ G.| C I h + ^ Ig ) 

dp t^ - G. tf 
■whom y ! — i 

“a h + '^1 d 

with Gp ™' Gp (1 +i 0 ) 

War])ing Rigidity of the Composite Beam 

Gp = E^ + Ep tp C^^ 

where C = warping constant of cross section 1 

as G„ in Appendix A . 
w 



where t 


( 



APPENDIX C 


DETAILS li'OR TflE FREE VIBRATION OF COMPOSITE 
BEAM IN COUPLED MODE 


C.1 Soction Properties 

The various section properties can he written 
for the composite cross section shown in Figure B-1 , as 


+ A. (t ^ - C.) 


-X) 1 c 


1-1 t^ _ _ 2 

— g — + A2 ( ^ 2 ^ 


I 


Ci 

lop = '<2 + h 

where I„ , \ C,,’ 's 

in previous appendices , 


C,2 Non Dimensional Terms in Equation (if.l8) 


2 

- g 1 

“2 = 

= g q 

“1+ - ® ^ 

where 

t 

8 = 

and Cg are the torsional and the warping 
rigidity given in Appendix B , 



c-2 


+ Pn A, 




1. I. + B_ 

) ^(2 Jj. i_£2) 


'1 


Pl^Ai + P^ A2 

¥horo Bp is the real part of . 


P., A., + P, A, 


c? ) X ( ^ 


E, + E^ I 


2 -"5, 


^2 


A^ + p^ Ar 


^0, + <’2 


+=52 15, 

) V C J- - ^ 

°1 h ^ ^2 


Noll climonsioiijxl. nabiu'al frequency is given by 


nj (1 + 1 V ) 


w (1 + 1 n ) 


1 ? 


+ E. 


1 


( p ^ A^ + p 2 Ag ) 1 


A 


vriicre w Is the natural frequency 

n in the nondimensional natural frequency 
c 

and n is the composite loss factor. 

V « n . 

G.3 Eloments of Matrix A in Equation (A. 20) 


'L'liu j; om of matrix A is same as given in 
A])ponc\lx A,>i'. itoro, ^ ^2 eighth root of 

cliarac tori otic polynomial equation given in Equation (A,13). 
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